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Abstract 

Retinitis pigmentosa (RP) is a group of inherited degenerative eye diseases char- 
acterized by mutations in the genetic structure of the photoreceptors that leads to 
the premature death of both rod and cone photoreceptors. Defects in particular 
genes encoding proteins that are involved in either the photoreceptor structure, 
phototransduction cascades, or visual cycle are expressed in the rods but ultimately 
affect both types of cells. RP is "typically" manifested by a steady death of rods 
followed by a period of stability in which cones survive initially and then inevitably 
die too. In some RP cases, rods and cones die off simultaneously or even cone death 
precedes rod death (reverse RP). The mechanisms and factors involved in the de- 
velopment of the different types of RP are not well understood nor have researchers 
been able to provide more than a limited number of short-term therapies. In this 
work we trace the progression of RP to complete blindness through each subtype via 
bifurcation theory. We show that the evolution of RP from one stage to another of- 
ten requires the failure of multiple components. Our results indicate that a delicate 
balance between the availability of nutrients and the rates of shedding and renewal 
of photoreceptors is needed at every stage of RP to halt its progression. This work 
provides a framework for future physiological investigations potentially leading to 
long-term targeted multi-facet interventions and therapies dependent on the partic- 
ular stage and subtype of RP under consideration. The results of this mathematical 
model may also give insight into the progression of many other degenerative eye dis- 
eases involving genetic mutations or secondary photoreceptor death and potential 
ways to circumvent these diseases. 

keywords: retinal degeneration, rod-cone RP, reverse RP, bifurcation set, 
RdCVF 
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Introduction 



Background of Retinitis Pigmentosa 

Retinitis Pigmentosa (RP) is a heterogeneous group of inherited disorders that affects 
approximately 1 in 4000 individuals (nearly 1.5 million) worldwide [IJ[2j[3]. This group 
of disorders is often classified as primary or syndromic RP. In the case of syndromic RP, 
which account for about 20-30% of all cases, the patients have an associated non-ocular 
disease with RP with abnormalities in one or more organs [U 12] • In contrast, primary 
RP involves only the eye. RP is primarily genetically programmed with more than 45 
genes associated with this disease identified to date. These genes account for only about 
60% of all cases. RP may be inherited from one or both parents (autosomal-recessive RP, 
autosomal-dominant, or X-linked RP, which account for approximately 20-30%, 15-20%, 
and 6-10% of the cases, respectively) or may be caused by a new mutation even though 
there may be no family history of the disease jU El El El [lj [7] . Less common modes of 
inherited RP include digenic and mitochondrial [2]. 

RP has been associated with mutations or errant genes encoding for proteins that 
are essential for rod survival and certain genetic products that are damaging to the rods 
[21 IH El E] . Even though mutations are normally only expressed in the rods their effects 
ultimately cascade to the cones. People usually are diagnosed only when they experience 
serious vision problems making it difficult to decipher the mechanism of development and 
the channels by which RP develops. Adding to this challenge is the variability of its 
progression and its time of onset within each individual [H [101 [11] . The sequences of 
events (i.e., clinical manifestations) of the various sub- types of RP are well-documented 
[21 H21 ELL] but this sequential progression varies from individual to individual. 

RP is often classified into various sub-types depending on the location of the vision 
field that is affected. The most common clinical manifestation is rod-cone RP and is char- 
acterized by difficulty seeing at night, followed by night blindness (nyctalopia), narrowing 
of the mid-peripheral field of vision (where the rods are more prevalent), and ultimately 
a loss of central vision [HI IH [3]. Night blindness is the hallmark of rod-cone RP but given 
that rod function is primarily affected, this is of no surprise since rods contribute solely to 
night vision. Even though rod-cone RP has a better overall prognosis and slower progres- 
sion in comparison with other sub-types of RP, the range of the onset of clinical rod-cone 
RP to its noticeable changes observed by the patient can vary widely in all cases. In this 
most common form of RP, the causal mutations affect genes exclusively expressed in the 
rod cells in an overwhelming number of cases and rod death far exceeds and precedes 
cone death [2] . Less common forms of RP see an equal loss of rods and cones or even a 
cone death exceeding rod death (cone-rod RP or reverse RP) [U12]- Cone-rod dystrophy 
(reverse RP) is associated with a cone loss preceding rod loss and as a consequence reduc- 
tion of the central vision, acuity, and color vision are manifested early in the progression 
of the disease, followed by reduction of the mid-peripheral visual field and night blindness 
at later stages. However, unlike the rod-cone RP where patients' color vision remains 
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good until the central vision is affected at levels of 20/40 or worse, patients with cone-rod 
RP experience deterioration of their color vision, ability to adapt to light and day vision 
very early in the disease as well as the loss of normal cone functioning, which leads to 
photophobia. 

Interactions of Photoreceptors 

Research that might shed light on improving the photoreceptors longevity has been 
paramount to physiological, anatomical, and medical research in the last 2 decades due to 
the prominent role that these cells play in the vision process [13] . The photoreceptors are 
the rods and cones, located in the retina, that convert light energy into nerve impulses 
through the process of visual transduction [H]. The rods and cones are specialized for dif- 
ferent aspects of vision: the rods and cones are different in shape (cylindrical vs. conical), 
contain different types of photopigment, and are distributed differently throughout the 
retina[15]. The reduction of normal light vision, visual acuity, detail perception, night vi- 
sion, and narrowing of the peripheral fields can be caused by a depletion of cones and rods, 
respectively [3]. The photoreceptors undergo tremendous oxidative stress over the course 
of a day and undergo a continuous cellular renewal and periodic shedding process of their 
outer segments (OS) in order to minimize depletion and help ensure proper functioning of 
these cells [T5l [T6| \T7\ [T8] . Photoreceptors shed approximately the same amount of mate- 
rial that is renewed over the course of a day so that the length roughly remains constant, 
with complete renewal in rods and cones occurring from 11 days up to 3 weeks depending 
on the photoreceptor and its location within the retina p3l [T9l [20| [21l [22| [23| [2^| [25| [26]. 
The shedding and renewal of photoreceptors, their interactions, and circadian rhythms 
were previously examined and discussed in greater detail [271 128] . 

Much research has examined the processes of renewal and shedding within the pho- 
toreceptors and suggests the photoreceptors together with the adjacent retinal pigment 
epithelium (RPE) operate as a functional unit PDIESII3DIQSIEIIE21E31IIEII311E31I221E51 
1361 12T1 1371 [38] . The rod and cone cells undergo a daily rhythmic shedding in which their 
photoreceptor outer segments (POS) are phagocytized by the RPE. The RPE recycles 
parts of the shed POS for future use by the photoreceptors and contributes additional 
trophic and growth factors utilized in the renewal process of these cells [291 El E3 EE] . The 
RPE transports, secretes, and redelivers essential growth factors such as PEDF, CNTF, 
bFGF, Gas6, PDGF, and ATP, metabolites (including glucose, retinal, and omega-3 fatty 
acids), ions, and water to and from the photoreceptors [381 SOI SH H21 EH 021 03]. It is 
unclear if other essential proteins are linked to the interactions of these two retinal tissues 
but the interdependency of the RPE and photoreceptors is key to the proper functionality 
of the eye. 

Pivotal to the vision is also the connection and interplay between the rods and cones. 
These cells have long been suspected of communicating directly or indirectly with each 
other [39]. It has been recently shown that the rods produce a protein referred to as 
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the Rod-derived Cone Viability Factor (RdCVF) that is essential for the functionality 
and survival of the cones [I5j [3j S6j H7J SHJ H9]. While the neural degeneration of the 
photoreceptors that typically results from genetic defects and mutations (such as RP) 
may disrupt and break the balanced shedding and renewal processes leading to shorter OS 
and, in some cases, the complete disappearance of the rods and cones [31 [50], the amount 
of RdCVF appeared to be the same in both normal and mutated rods prior to their 
degeneration. Additionally, a subsequent decrease in RdCVF was reported to be strongly 
correlated with the decrease in the number of rods during their degenerative process [4"6] . 
The effect of RdCVF appears to be independent of the casual mutation leading to rod 
degeneration [IHl HHJ HE] . Experiments have revealed that RdCVF had a significant cone 
rescuing effect, that its protective effects did not extend to degenerating rods, and that 
the presence of the rod cells is not necessary for cone survival as long as RdCVF is present 
[4"8| |4"5[ |4"6| Wf\ |4*9[ [28] . Moreover, RdCVF was also shown to significantly preserve cone 
function, a very promising discovery since maintaining functional cones even when 95% of 
the cones are gone may prevent blindness [IB]. The discovery that rods produce RdCVF, 
which is essential for cone functionality and viability, and that their death deprives the 
cones from this element and triggers their degeneration has far-reaching implications. It 
is yet unclear whether the RPE plays any role in the transport of this protein but it is 
clear that RdCVF provides an essential direct help from the rods to the cones without 
negatively affecting the rods. 

Hypotheses of Photoreceptor Degeneration and Death in RP 

It is crucial to identify the precise cause(s) of both rod and cone cell death, as it is 
the secondary wave of photoreceptor death in RP patients that results in total blindness 
[29[ IBTj . Substantial experimental research (consisting mainly of animal models) has 
been devoted to this effort. There are several hypotheses/theories as to what leads to 
photoreceptor death in RP, but none of them have been definitively proven [29]. All 
theories suggest indirectly or directly that mutations result in gene products that are 
damaging to photoreceptors (including their signaling pathways and cell-cell interactions) 
or that mutations cause lack of a protein that is needed for photoreceptor survival [Hj. 
Among these theories is one which led to the development of antioxidant therapy for cone 
survival and is based on the hypothesis that the death of rods leads to higher oxidative 
stress in the retina, which eventually kills the cones [29]. This was demonstrated by 
showing that antioxidants promote cone survival and vision in rdl mice as they get older. 
Another theory comes from the idea that the rods produce some kind of signal that is 
essential for maintaining cone viability, so that the disappearance of rods for whatever 
reason would deprive the cones of this signal and trigger their degeneration [3]. This theory 
was supported through various experiments on mice and rats and led to the discovery of 
the RdCVF. A third theory comes from the synergistic interactions of the RPE and the 
photoreceptors since defects in genes normally expressed in the photoreceptors can lead 
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to degeneration of the RPE and vice versa [38]. The inability of the RPE to secrete 
neuroprotective and growth factors such as PEDF, bFGF, and CNTF that can protect 
photoreceptors from damage has been implicated in the pathogenesis of retinal diseases 
( [38] and references therein). Inappropriate signaling between photoreceptors and bipolar 
cells whose signaling pathway remains intact for RP patients has also been postulated 
as a potential reason for the secondary death wave of photoreceptor in rod-cone RP 
patients [52]. Toxicity released by mutated and degenerate rods is another hypothesis 
but this theory has been challenged by experiments which showed that cone death is 
uniform in space irrespective of where they are in relationship to the mutated rods [53| 13]. 
Another related hypothesis is that microglia activated by ongoing rod death migrate to 
the photoreceptor layer and secrete toxic substances that kill cones [5"3] . 

Depletion of photoreceptors or abnormally short OS that result in oxygen toxicity only 
explain part of the cause of cell death in RP [551 [13]. Although rods die by apoptosis, which 
generally avoids release of potentially toxic cellular metabolites and prevents extension of 
cellular breakdown to surrounding regions of tissues, it seems possible that degeneration 
of the abundant rods (numbering about 90-120 million in the human eye) could have 
adverse effects upon the adjacent cones and nearby bipolar cells [3]. Even though there 
are some modifications in the inner retinal neurons these nearby neurons do not undergo 
a widespread death following the first wave of photoreceptors death, as one would expect. 
The main flaw in these toxicity theories is that cone cell death occurs very slowly and 
continues long after all rods have died. While toxicity might explain photoreceptors 
degeneration in the short term (at the early stage of the disease) it fails to hold in the 
long run [3]. 

The theory of oxidative stress fails to completely explain the sequential photoreceptor 
death in reverse RP (or in cases when both cones and rods die simultaneously) and in 
the rod-cone RP cases it raises the questions as to why rods die preferentially and why 
only cones are saved when treated with antioxidants. In an attempt to address this 
researchers have postulated that rods might be more sensitive to reactive oxygen species 
(ROS) and thus they die first, which then probably further increases the oxidative stress 
in the retina. The generation of ROS, by proteins in various retinal layers, is part of 
the normal maintenance of the retina [29J. Rods are more numerous than cones and are 
metabolically active cells with a high level of oxygen consumption. Choroidal vessels are 
not subject to autoregulation by tissue oxygen levels and as death of rods occurs, the 
tissue level of oxygen in the retina increases [561 EZ]- Oxygen levels in the outer retina 
vary inversely with oxygen consumption; as the number of rods in a particular area of the 
retina decreases, oxygen consumption decreases, and thus oxygen levels increase. Based 
on this, one would expect regions of human retinas with the highest rod/cone ratio (the 
far periphery of the eye) to be the first regions to experience photoreceptor death due to 
their exposure to dangerously high levels of oxygen [8] . This is what we see in the case 
of rod-cone RP where the peripheral visual field is constricted first but not in the case 
of reverse RP. While this theory might not depict the entire picture it provides one of 
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potentially multiple causes of cell death in RP as it has been established experimentally 
that high tissue levels of oxygen are damaging to photoreceptors and result in cell death 

PIES]. 

Another valuable insight for the death of rods and cones in RP comes from our current 
understanding of the RPE as a vital component of the visual system that supports normal 
photoreceptor function. Integrated interactions between photoreceptors and the RPE are 
key in the visual process. Among these interactions are the visual cycle (in which the 
structure of retinal is changed by light transduction: in photoreceptors from 11-cis- to 
all- trans-retinal and in the RPE from all-trans- to 11-cis- retinal) and the phagocytosis of 
photoreceptor OS (which contributes to the photoreceptors shedding process and the op- 
timal photoreceptor structure necessary for light perception). Experiments have revealed 
that signs of photoreceptor apoptosis due to intensified light exposure were fully reversible 
potentially due to increased secretion and synthesis of different neuro protective factors 
by the RPE such as PEDF, CNTF, and bFGF, which are known to be protective against 
light damage Among the essential proteins linking the RPE and photoreceptors is 
the interphotoreceptor retinoid binding protein (IRBP) which is primarily a product of 
the rods, with less of a contribution from the cones [601 |6T| |62| [3j. A decrease in IRBP 
due to the death of rods as well as a loss of the RPE-recycled products of rod OS shed- 
ding may contribute to the observed cone OS shortening and ultimately, the cone cell 
death found in RP retinas. This is mainly attributed to an associated local deficiency in 
vitamin A, as it is known that vitamin A deficiency causes photoreceptor OS shortening 
and cell death [63j [M]. The inability of the RPE to provide essential neuro protective 
and growth factors to individuals with RP opens a realm of possibilities that call for a 
thorough investigation of the interactions of photoreceptors and their exchange of trophic 
factors. 

Another protein that seems to be key in the interactions of the photoreceptors respon- 
sible for the viability of the cones is the RdCVF. RdCVF is a truncated thioredoxin-like 
protein specifically expressed by the rods. Experiments demonstrated a 40% rescue ef- 
fect in cones and thus slowed down the degeneration of cones in chick and mouse models 
[45"| HB] . Additional experiments revealed that RdCVF had a bigger effect in the function- 
ality of cones than in their viability [35]. Although RdCVF is a good candidate to target 
for cone survival in RP, a number of issues remain to be addressed. For example, the 
mechanism of action underlying the cone survival effect needs to be better understood as 
well as the intracellular signaling pathways that are involved jH]. Another issue that is 
currently being researched is the delivery mechanism of RdCVF as injection of a protein 
but this might not to be a satisfactory option for the long term. Gene-based delivery 
through viral vectors or using encapsulated cell technology would ensure a more stable 
and persistent release of RdCVF [J5]. Irrespective of the mutation, the final common 
pathway of the first wave of photoreceptor death is degeneration through apoptosis; the 
mechanism underlying secondary wave of cone (in rod-cone RP) or rod (reverse RP) death 
remains unknown [3J. 
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Mathematical Model 



In an attempt to shed light into the underlying mechanisms and factors involved in the 
development of the different types of RP we create a mathematical model of photoreceptor 
interactions that explicitly incorporates RdCVF and implicitly includes other essential 
photoreceptor neuroprotective and growth factors secreted by the RPE. We consider two 
different classes of rod cells, mutated and normal, to represent those that have had their 
functionality affected by the genetic mutation(s) and those that have not [2H1 ESI ES] (see 
Figured] for a simple physiological illustration). This distinction between the two classes 
of rods is essential due to the numerous genes and processes affected by the mutations 
associated with RP[12] [2] [8] and the great variability of the disease [67]. While much 
RP research has characterized the different types of mutations that cause RP, we are 
not modeling and do not address how the mutation occurred or what type of mutation 
it is or what it affects. Rather, our model focuses on the interaction and structure of 
the photoreceptors and their OS in relation to the renewal and shedding process. We 
simply assume that there is a mutation that has already been expressed and therefore 
some of the rods cells can be thought of as "mutated rods" [67]. We similarly do not 
model the visual transduction or retinal metabolism [13]. Thus our model can be thought 
of as representative of the photoreceptor and RPE interactions at or after the onset of 
RP[2B]- As is often done, we used previous mathematical models of the eye together with 
experimental evidence to assume that both eyes are identical and thus consider only one 
eye [651 1211 EH1 12] • This work builds on our previous mathematical model of photoreceptor 
interactions in a healthy eye but now we incorporate mutated rods that are present in RP. 
For the sake of mathematical simplicity, we also ignore any spatial effects, for example, the 
different spatial distribution of the photoreceptors at various locations in the retina [66] . 
Two other mechanistic models of RP focus on the kinetics of neuron degeneration and 
lateral inhibition of nearby cells and not on cone and rod photoreceptor interaction [66| 169] . 

In order to focus on the mechanisms behind the various manifestations of RP, we 
create a mathematical model that includes two different classes of rod cells, R m and R n 
(corresponding to mutated and normal rod OS, respectively), one type of cone OS (green, 
red, blue cone OS are grouped together), C, and a trophic pool (representative of the 
RPE's neuroprotective and growth factors) from which the rods and cones draw their 
nutrients, T. We utilize the framework from epidemiology and population dynamics to 
model the interactions between the photoreceptors and trophic factors as a predator-prey 
system; see Figure [TJ Even though no new rods or cones are created in a mature retina, 
the process of shedding and renewal lends itself to the interpretation of birth and death of 
the rods and cones and we incorporate this interpretation into our model [32], [70], [33]. We 
view our units as "cells" and therefore fractions of a cell are the POS in these three types of 
photoreceptors. Experiments have shown that in a rod-free retina the absence of cone OS 
corresponds to blindness even though the living cone cell remains [7T| 126"]. We will take this 
same perspective. Similarly the nutrients represent fractions of the trophic pool consisting 
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of growth factors, metabolites, ions, and water (e.g., PEDF, CNTF, bFGF, Gas6, PDGF, 
ATP, glucose, retinal, and omega-3 fatty acids) to and from the photoreceptors provided 
or secreted by the RPE [38l ffifl HU 021 E3 SSI IB] . In RP, degeneration can also extend 
to the outer retina nuclear layers and the ganglion cell layer but we focus this work only 
on the photoreceptor cells [ZD U2\- Our model extends the healthy eye model of [28] : 

R n = Rn(a n T - fi n - m) 

R m = Rm(a m T — fi m ) + mR n 

C = C (a c T — fi c + d n R n + d m R m ) 

f = T(T — kT — f3 n R n — f3 m R m — 7C), (1) 

where all parameters are nonnegative; see Table [2] for definitions of the parameters. We do 
not distinguish between the various factors that constitute the trophic factors so as to focus 
on the critical interactions of the rods and the cones [28J. Furthermore, this allows our work 
to stay consistent with all experiments in the literature where no distinction between a rod- 
trophic and a cone-trophic pool is made [151 129] . Our model assumes that both the rods 
and the cones take more nutrients from the trophic pool than they contribute (as recycled 
products) and that they convert these trophic factors to new OS discs at some ratio. To 
the best of our knowledge, no experimental work has been done to date that discusses 
the production rate of trophic factors in the absence of photoreceptors; this current work 
supposes that the production of the trophic factors is logistic in nature and thus will reach 
an equilibrium (governed by T and k) in the absence of photoreceptors, instead of growing 
without bound [2H]- The structure of this mathematical model incorporates the presence, 
through the parameters d m and d n , of the RdCVF protein that helps cone survivability 
and functionality yet does not negatively influence the rod population [35l [3l [H 08]. 
Experimental evidence suggests that the RdCVF directly affects the cones and thus we 
model this as a direct rod-cone interaction, the necessity of which was already established 
in [28] . The literature does not suggest that mutated rods stop or decrease production of 
RdCVF and thus our model includes this possible contribution of RdCVF by the mutated 
rods via d m . This model reduces to the previous work in [28] if we ignore the terms related 
to the mutated rods, i.e., m = R m = d m — (3 m = 0, and do not consider saturation of the 
trophic pool, k — 0. 

We define the following terms: 

D m = ^, D n= >^, D c= '^ D T= T -, (2) 

which are key in our analysis. Based on the equations ([T]), we determine steady state 
solutions (R n , R m ,C,T) that correspond to the various physiological stages of RP and 
whose corresponding stage and subtype of RP is summarized in Table El Our system 
yields seven equilibria given by 
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(3) 



Ax = a n a m (D m - D n ), A 2 = a n a c (D c - D n ), 

A 3 = a n k(D T - D n ), Bx = a m k (D T - D m ), 

B 2 = a c k(D T - D c ), B 3 = a m a c (D c - D m ). 

The ratios defined by the Di for i — m,n,c,T represent key quantities that govern the 
progression of RP and will be discussed in Appendix 1. 



Progression of RP through Equilibrium Solutions 

Analysis of the seven steady state solutions Q give conditions on the parameters for the 
linear stability of the various stages of RP [39] . Once a particular stage of RP has been 
reached, these conditions represent quantities that must persist in an individual to prevent 
further progression of RP (e.g., moving to a more advanced stage of RP) and equivalently 
they give the parameters involved if RP progression continues. These conditions involve 
the terms in <^fy: D m , D n , D c , and Dt, which turn out to be the key quantities predicted 
by our model to govern the progression ofRPQ The first three of these terms are the ratios 
of the OS departure (due to shedding and/or mutation) to the renewal for a given class 

1 A non-dimensionalization of the original system, while reducing the parameters from 14 to 10 yields 
the same 4 essential ratios in © that are involved in the various bifurcation sequences, without additional 
insight gained to date with the non-dimensional system. Thus, in order to keep the problem tied to the 
physiology as much as possible, we leave our system in its original dimensional form and focus on the 
relationship between the Di as the mechanisms governing the progression of RP. 
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of photoreceptors. The fourth term, Dt, is the ratio of the natural growth rate to the 
carrying capacity of the trophic factors in the absence of any photoreceptors (complete 
blindness). Our results are supported by previous experimental research suggesting that 
disturbances or imbalances of disc shedding to disc formation is associated with RP [73] . 

Linear stability analysis of the steady state solutions E^ % — 1, • • • , 7 allow us to de- 
scribe the progression of RP for various subtypes within rod-cone RP and reverse RP. 

Result: We have observed the equilibria of system (QJ) in multiple trans critical bifurcations 
yielding the following progressions: 

1. E 7 -> E A -+ E 2 ; E 7 -> E 5 -> £ 4 -> E 2 (rod-cone RP) 

2. E 7 ^ E 6 ^ E 2 ; E 7 ^ E 6 ^ E 3 ^ E 2 (cone-rod RP) 

3. E 7 — > E 5 — >■ E 3 — > E 2 (cone-rod RP or simultaneous photoreceptor death) 

4- E 7 — >■ E 2 ; (simultaneous photoreceptor death) 

In each of 1-4, a given equilibria Ei is stable and loses its stability through a transcritical 
bifurcation (—$■) as the subsequent solution becomes stable. The stability of the equilibria 
depends on D m , D n , D c , Dt- 

While not all of the transcritical bifurcations are obvious from inspection, calculations 
in a computer software package such as Maple show that at each bifurcation, the two 
equilibria become one as they switch (at least) one direction of stability for realistic 
parameter values. When D n = D m = D c = Dt in parameter space, equilibria E 2 , ■ ■ ■ ,E 7 
coalesce in an equilibrium point with three zero eigenvalues. We have examined the flow 
in the corresponding three-dimensional center manifold of the non-dimensional version 
of dl]) but have not gained any additional insight regarding the observed bifurcations. 
The numerical results on the center manifold flow have confirmed the numerical results 
in the original system regarding the stability of only one equilibrium being relevant and 
stable at any given point in the realistic range in parameter space. However, we have 
observed stable periodic solutions for certain (unrealistic) parameter values. Following 
these solutions with MatCont has not yet allowed us to conclude the bifurcation that 
gave rise to these limit cycles (one exists in the full 4-dimensional space while the other 
exists in the 3-dimensional subspace given by R n = 0). This aspect of the work is still 
ongoing. Thus, while we have analytically observed the sequence of bifurcations in the 
above Result for specific realistic parameter values, we are only able to conjecture that 
this result holds for all choices of parameter values. The statements below are thus certain 
for realistic parameter values and conjectured to hold in general. See Appendix 1 for a 
discussion of the analytical results regarding stability of the equilibria, Appendix 2 for 
a discussion about parameter estimation of the model, and Appendix 3 for other stable 
solutions observed. 
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It is important to observe that E4, for example, represents a state in which both types 
of rods have died off and only cones remain; however, the number of cones present in a 
given E4 value may differ widely and depends on the parameters. Thus, the bifurcation 
E 7 — > E4 allows for the possibility that the cone levels in E 7 have remained the same 
as the rods die off (as the values in E 7 get closer to those in E4 as the parameters are 
varied), but also allows for the possibility that some/many of the cones have also begun 
to die off. Thus while we can think of E4 as the typical state in rod-cone RP in which 
all rods have died off while the cones remain, we must remember that the level of cones 
in E4 may not be the original numbers of cones seen in a healthy eye. Keeping this in 
mind, we can interpret the sequences of transcritical bifurcations of the Result as various 
subtypes of RP by asking "which photoreceptors remain just after the bifurcation" as we 
vary a given parameter or set of parameters. 

The sequences given in 1 describe the progression of rod-cone RP via two different 
paths. The transitions E 7 — > E 4 — ?► E 2 might represent a more aggressive subtype of RP 
in which night blindness occurs because of a simultaneous loss of mutated and normal 
rods (E 7 — > E4) before the complete death of cones (E4 — > E 2 ). The sequence E 7 — > E 5 — >■ 
E4 — > E 2 might be a less aggressive subtype as the normal rods die off first (E 7 — > E5) 
followed by the mutated rods (£5 — > E4) before the cones die off. In both sequences, we 
have an intermediate state in which only cones remain before eventually dying off. The 
sequences given in 2 describe the progression of cone-rod RP via two different paths, with 
the latter having the interpretation of both cone-rod RP and simultaneous death of rods 
and cones depending on the actual equilibrium values. The transition E 7 — > E e — > E 2 
might represent a more aggressive form in which cones die off before either normal or 
mutated rods die off (E 7 — > E 6 ) and then both types of rods die off together (E 6 — > E 2 ). 
In the path E 7 — > E 6 — ► E 3 — > E 2 , the cones are predicted to die off first [E 7 — > E 6 ), 
then the normal rods (E$ — > E 3 ), followed by the death of the mutated rods rendering 
complete blindness (E 3 — > E 2 ). The sequence E 7 — > E 5 — > E 3 — > E 2 (in 3) might 
represent a somewhat less aggressive form of cone-rod or simultaneous death RP in which 
only normal rods die off (E 7 — > E 5 ) before the cones completely disappear (E 5 — > E 3 ). In 
this sequence the individual will have limited night vision and thus clinically this would 
correspond to a cone-rod RP as long as night vision persists to some degree. In these 
three sequences (of 2 and 3), it is a type of rod that is the last photoreceptor to die off 
thus leading to the classification of all three as cone-rod RP. In contrast, the sequence 
given in 4, E 7 — > E 2 , describes an additional sub- type of RP in which cones and rods 
die off simultaneously. If the functionality of the rods has been extremely affected by 
the mutations to the point where night vision has been completely impaired such that 
E 3 contains no functional rods at all then the sequence E 7 — > £5 — > E 3 — >■ E 2 can be 
considered to be in category 3 where clinically cones and rods are considered to die off 
simultaneously. For the manner in which RP is addressed in this work photoreceptors 
that contained OS that shed and renew are considered to be functional to some degree. 

It is important to again note that some of the above sequences, for example, E 7 — > 
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Ei — > E 2 may in some situations be clinically classify as the same subtype as E 7 — > 
E 2 if the cone levels of the E± equilibrium are very close to zero and disappear very 
shortly after the rods disappear. To the best of our knowledge, no experiments have 
made the distinction between normal and mutated rod death during the progression of 
any sub-type of RP nor has any research been able to follow the complete progression 
of death for a given type of photoreceptor (i.e., produce a time series for the cones or 
rods in a given mammalian subject). The results predicted by this model thus provide 
a potentially crucial understanding of the progression of RP in terms of the parameters 
that give rise to these stage changes. This knowledge gives the framework that will guide 
future experiments focusing on potential therapies to slow the inevitable progression of 
RP. Insights into what must be considered comes from the stability conditions of the 
transcritical bifurcations. 

Since the terms D n , D m , D c , are all involved in the conditions for the stability 
and/or physiological existence of equilibria E 2 , ■ ■ ■ ,E 7 , future experiments exploring po- 
tential long term therapies and interventions would need to keep in mind these quantities 
and their interrelatedness as these yield the multi-faceted therapies that can halt RP 
brought about by the change in these quantities. By varying between one and three pa- 
rameters of the system ([I]), we can thus observe many different progressions of cone- rod 
RP and rod-cone RP (see Figures [2] and [3]) . Our model suggests examining experimen- 
tally the causes of the shedding and renewal imbalances and of changes in the trophic 
pool, mathematically seen as changes in D m , D n , D c , and Dt, and studying how to alter 
these quantities in RP patients [73]. In our system, we have observed exactly one stable 
equilibrium Ei, % = 2, • • • ,7 that is stable for a given choice of realistic parameters. This 
limits the possibilities for individuals with RP and thus it becomes even more crucial to 
have a better understanding of how to clinically control the ratios D m , D n , D c , and D T . 

Conclusion 

The work presented here offers the first quantitative model of photoreceptor interactions 
that can track the progression of RP through all the various stages and subtypes of this 
disease and thus offers alternatives to stop its evolution at each stage. Our model predicts 
the progression of rod-cone RP (E 7 — » £4 — > E 2 and E 7 — > E 5 — > E± — )■ E 2 ) and cone-rod 
RP (E 7 -> E 6 E 2 , E 7 ->■ E 6 — )► E 3 ->■ E 2 , and E 7 ->■ E 5 -)■ E 3 ->■ E 2 ) by examining 
the ratios in the quantities D m , D n , D c , and Dt- Because the stability of the equilibria 
are governed by these ratios, our model suggests that future physiological experiments 
and research should focus on the shedding and renewal for each class of photoreceptors 
as well as the carrying capacity of the trophic factors as potential leads into long-term 
treatments for RP [73] . 

The stability of the various stages of RP depend on the difference of these ratios, 
Di — Dj, for i ^ j and i,j = m, n, c, T given above. The steady state solutions, represen- 
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tative of the RP stages, become physiologically relevant through transcritical bifurcations, 
when two steady states coincide and then switch their stability [71]. Analytical and nu- 
merical evidence indicates that only one of the steady state solutions is stable for realistic 
parameter values; for certain unrealistic parameter values, bi-stability has been observed 
(see Appendix 3). By varying two parameters, we can track the progression of RP as we 
cross the regions in which a given steady state solution representing the current stage of 
RP becomes unstable; see Figure [2j 

Progression of RP culminating in complete blindness means we begin in the region 
where Ej is the stable solution (first stage of RP) and end up in the region where E2 is 
the stable solution (complete blindness). Sample paths 1 and 2 in Figure [2] illustrate the 
variability in RP and how identical changes in the parameters in two hypothetical different 
patients (both with realistic parameter values) could yield different clinical manifestations 
and even sub-types of RP[TTl [8j [10]. In path 1 we track the progression by varying the 
same D m , D n , and values of two hypothetical individuals with all other parameters 
fixed at the same values except different D c values (2.34 x 10 6 and 2.45 x 10 6 , respectively 
with both realistic) and observed very different manifestations of RP, Ej — > Eq — > E2 
versus E 7 — > E 2 . In sample path 2 we see that while the observed sequence of RP stages 
might be the same, the range of parameter values in each RP stage may vary drastically in 
two patients [TBI [76] even when only one parameter is different in their visual system. The 
different ranges are indicative of the different options of available to halt the disease. This 
finding is consistent with the established variability of RP between individuals [8j [761 143] . 
Furthermore, this result opens the opportunity for multiple individual targeted therapies 
to slow the progression of RP. 

While we may be able to see all possible stages of RP in an appropriate parameter 
region by varying at most two parameters, it is unlikely that this will be the norm phys- 
iologically since these biological components are all interconnected and disturbances in 
one can trickle down to others [TO]. Thus Figure [3] shows how we could progress from 
Ef to E4 to E2 in the most common rod-cone RP by varying three parameters, thereby 
varying D m , D c , and D?- The relative size of the stability regions of E 7 and E4 in this 
space illustrate the various ways to remain in the RP stages corresponding to E? or E± 
both of which permit day/color vision. The progression of rod-cone RP illustrated in the 
Sample Path (Figure [3]) requires a decrease in the ratio of shedding to renewal of the 
cones (D c ) together with a decrease in the necessary trophic factors (i.e., nutrient star- 
vation) for the vision system of this hypothetical individual to fail to meet the metabolic 
and functional needs of the photoreceptors (in the particular the cones). The changes in 
shedding, renewal, and trophic factors required to move from one stage to another depend 
on the differences of the relevant quantities D{ — Dj, i ^ j. Counteracting these currently 
inevitable trajectories of rod-cone RP in our model translate to changing D c and such 
changes correspond to new therapeutic approaches that bring back a balance between 
shedding and renewal of the cone OS [73] . 

Currently there is no standard long-term treatment for RP, in part because the un- 
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derstanding of RP is incomplete. There are many different theories for the physiological 
processes triggering photoreceptor death and RP progression. Some of these causes, while 
still having physiological value, can be contradictory or can explain only certain aspects of 
RP progression, and/or only hold for the early stages of the disease. Research in this area 
has improved the understanding of RP tremendously yet the possible mechanisms that 
trigger its inevitable progression and offer potential long term therapies are still open ques- 
tions. Our model provides a possible unified understanding to these two open questions. 
Recent experimental research has focused on nutritional and drug therapy, the visual 
transduction cascade, and the renewal and shedding of the POS [2[ [TTJ [39], H31 HEJ 1751 138] . 
For example, high dosage of vitamin A supplements have been utilized as a nutritional 
therapy as it has been shown to cause an increase in the trophic factors of nutrients es- 
sential for proper photoreceptor functioning [75] and our model illustrates this and similar 
nutritional therapy in addition to growth factor delivery by an increase in Dt- Likewise, 
changes in the shedding and/or renewal of the POS, which have been proven to be critical 
in the maintenance of healthy photoreceptors and are part of our main results, correspond 
to changes in D m , D n , or D c in our model depending on the specific type of photorecep- 
tor that is affected. Gene therapy, which is currently being explored experimentally, in 
our model corresponds to an impact on D n directly by replacing mutant or errant genes 
that affect normal rod shedding and/or renewal, and perhaps the others indirectly if their 
renewal or shedding rates were affected through changes in D n [4T] . 

A promising new experimental approach to halt RP (and prevent total blindness) 
corresponds to maintaining the stability of E 4 (when only cones remain) in our model 
and it involves the discovery of the RdCVF [I5l [91 [391 SSI H7J- Mathematically and phys- 
iologically, RdCVF is necessary in order for the cones' viability and thus injections of 
RdCVF offer a rescue effect on cone survivability and functionality in the absence of 
rodsjlH]. However, other approaches to avoid complete blindness arise from mathemati- 
cally analyzing the stability of E± and the components that make it stable, which yield 
the conditions D c < D n , Dt > D c , and D c < D m . While D c is common in each and 
could be addressed by changes in the cone shedding or renewal, the other three quantities 
also must be considered to ensure that the cones-only steady state, E4, remains viable. 
A myriad of experimental treatments could be investigated to address these conditions. 

RP is the term used to describe a collection of inherited degenerative eye diseases, 
which are not fully understood and for which reliable long term treatments still await j2[ 
[T2[ [3j [HI [791 HQ]. While the vast majority of mutations have been identified in the genes 
that affect the rods, both rods and cones ultimately are affected [5[ IT[ IT2 [ I3"9" [ I75 [ ITT] . There 
are numerous other degenerative eye diseases that are not termed RP but that follow 
similar mechanisms of cell death triggered by mutations. Given the general structure 
of our model examining photoreceptor interactions, we believe that this work will help 
spur mathematical and experimental research in this area with the analogous conclusions 
inferring multi-faceted approaches to slow the progression of many other degenerative 
diseases. 
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Appendix 1: Stability of Equilibrium Solutions 

In order to examine the stability of the steady state solutions (J3J), we consider the linear 
stability of each, together with the numerical observations (e.g., see Figures [2]l3]) : 

Ex = (0,0,0,0) 

We can easily calculate the eigenvalues Aj = — (fi n + m), —fi m , —fi c , I\ which shows that 
Ei is always a saddle point (unstable). 

E 2 = (0,0,0, D T ) 

This equilibrium is always physiologically relevant and denotes the trivial case in which 
there are no rods (of either type) and no cones but only trophic factors. Physiologically 
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this corresponds to total blindness in an individual. The eigenvalues are again readily 
calculated as Aj = \B\, \B 2 , —T. Any of the sub-types of RP could eventually lead 
to this state. 



E 3 = (0 



0,D m ) 



The equilibrium E 3 represents the long term existence of mutated rods and trophic factors 
in the absence of cones and normal rods. E3 is physiologically relevant when B\ > 0, 
which is exactly one of the conditions for E 2 to be unstable. Thus as E 3 ceases to 
be physiologically relevant, it loses stability as E 2 gains its stability in a transcritical 
bifurcation, illustrated by the final stages of the sequences E 7 — >• E e —> E 3 — >■ E 2 and 
Ej — > E5 — > -E3 — > E 2 . The eigenvalues of £3 are given by 



This will be stable for A±, C| < and i?i > 0, the latter also being the condition for E3 
to be physiologically relevant. 



E4 represents the situation where an individual only has day vision. It illustrates the 
long term existence of the cones and trophic factors with both normal and mutated rods 
having disappeared. It is a physiologically relevant solution when B 2 > 0, which is exactly 
one of the conditions for E 2 to be unstable. The eigenvalues are given by 



This equilibrium point corresponds to a stage in the typical rod-cone RP in which both 
types of rods die off before we observe the secondary wave of cone deaths (months and 
even years after all rods have disappeared) [8] . Thus, we would expect this equilibrium 
point to be a temporary state in an RP patient that physiologically ceases to exist when 
the natural carrying capacity, k, of the trophic pool falls below a threshold value. This 
will render the pool unable to sustain the demands required to maintain cone life and 
functionality (i.e., nutrient starvation) when B 2 becomes negative. The functionality of 
the photoreceptor in this model is illustrated by a balance between natural renewal, at, 
and shedding, //j,z = n,m,c of the OS. Physiological experimental results have illus- 
trated a lack of functionality by shorter OS and a larger diameter of the disc at the tip 
of the OS Disturbances in the shedding or renewal process can lead to these mor- 
phological changes. As a result of disturbances on this delicate balance the cones would 
also eventually die off [73]. The stability of this equilibrium, E4, may give insight into 
why the secondary wave of cone death may follow. Physiologically, we note that, e.g., 
D r = — can increase with an increase in the shedding rate or a slowing of the renewal 




E 4 =(0,0,^-,D c ) 
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rate. As D n , D m , D c , Dt are all involved in the stability and/or existence of this equilib- 
rium, experiments would need to keep in mind the potential interplay of these quantities 
in determining exactly what processes may effect slight but crucial differences in these 
quantities. 

The eigenvalues of E 5 , E e , E 7 are too cumbersome to write explicitly but we can still 
gain insight by examining where each is physiologically relevant and in relation to the 
stability of the other equilibria through bifurcations. We can, however, write the explicit 
formulation for any A = bifurcations. Because our Result states that the bifurcations 
describing the progression of RP occur through transcritical bifurcations, these expres- 
sions will be useful. 



E 5 represents the long term existence of mutated rods, cones, and trophic factors with only 
the normal rods having disappeared. This point is physiologically relevant when > 
and i?3 > 0. Remaining at this state is desirable for someone with RP, as presumably 
the patient will still have color and day vision (due to the presence of the cones) and 



Equilibrium point E 5 corresponds to an earlier stage of the typical rod-cone RP in which 
only the normal rods have died off and both cones and mutated rods still remain. We 
would again expect this equilibrium point to physiologically represent a temporary state 
before cones and mutated rods eventually die off as parameters in the system change. 



This equilibrium represents the death of all cones with the normal and mutated rods as 
well as the trophic factors remaining. Besides depending on the renewal and shedding 
rates, E$ also depends on the rate at which both types of rods remove trophic factors 
from the trophic pool (/3 m ,/3„) and the production of the trophic factors itself (r). This 
equilibrium can be interpreted as the first stage in the progression of reverse RP in which 
an individual loses color and day vision. The eigenvalues are given by 

C^(d n Ai + d m a n m) 
m(3 m a n + (3 n A 1 

and a complicated cubic polynomial. Any A = bifurcations are given by 

(/i n + m) 1 A 1 A 3 C; _ 





E& — (Rle, RmQ, 0, D n ) 
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In the paths involved in the various progression of RP, E§ loses stability in a transcritical 
bifurcation with E 3 in the sequence E 7 — > E$ — > E 3 — > E 2 or with E 2 in the sequence 
E-j E§ E 2 . Similarly, Eq gains stability through a transcritical bifurcation with E-j. 
Understanding the mechanisms (i.e., the parameters involved) in Ej — > E§ is key as this 
would give insight into circumventing any vision loss. 

This equilibrium represents the only state in which all photoreceptors continue to survive. 
It illustrates the first stage in RP where potentially there are no clinical manifestations of 
RP. E 7 , in contrast to E e , besides depending on parameters governing production of the 
trophic factors additionally depends on the ability (or lack thereof) of the cones to utilize 
the RdCVF, quantified through d n and d m . The eigenvalues are given by a complicated 
4th order polynomial. The A = bifurcations are given by 

(fj, n + m^A^C; 

This is the most desirable state for an individual with RP — both types of rods, the cones, 
and the trophic factors are all still present in the system. Through an in depth analysis 
of this steady state, we can examine how to prevent RP from running its normal course 
by considering the conditions that will bring about qualitative changes that will render 
this state unstable and move the system to another steady state previously mentioned. 



Appendix 2: Parameter Estimation 

In obtaining estimates on the parameters for (CQ), we begin with the healthy eye (non- 
disease) version 

Rn = Rn(d n T — fl n ) 

C = C(a c T - ji c + d n R n ) 

f = T(T — kT — (3 n R n — 7C). (4) 

A few key experiments allow us to estimate five of these parameters. We then solve the 
inverse problem with 9 data points from the aforementioned experiments to estimate the 
remaining four parameters. 

Shedding of photoreceptor outer segments 

Numerous experiments, some in vivo and some in post-mortem patients, have estimated 
the number of rod OS and cone OS that are shed per day. This number is typically given 
in the number of days necessary to regenerate, i.e., all the OS to shed and be replaced by 
new ones. Estimates depend on both the location of photoreceptor and the species under 
consideration. If we consider only rhesus monkey and humans, we obtain an approximate 
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shedding range (in days) of fi n E [|, ^] ,/i c 6 [g^, ^] . See [15j [22J [THJ [80] for some 
of these key experiments and results. 
Rod-derived Cone Viability Factor 

The ground breaking experiment of Leveillard et al. gives insight the approximate value 
for the help the rods give the cones [45J. In their experiment in the rdl mouse they showed 
that the presence of RdCVF in a nutrient rich medium resulted in a 40% rescue effect. 
The number of cones per mm 2 after one week was 4200 compared with 3200 in their 
cultured retinal explants. Assuming an exponential decay to represent the degradation in 
the absence of RdCVF, we consider the equations 

C = -bC, C = -bC + d n R n C (5) 

to represent the cones in the absence and presence, respectively, of RdCVF. Solving and 
eliminating b and assuming that we have a level of RdCVF remaining near constant levels 
(where we let _R max represent the number of rods that would be needed to generate this) 
yields 

lnC^) 

d n R™ x « ^ff (6) 

Approximating _R max as the average number of rods per mm 2 (87.0 x 10 3 [T5l 181] ) gives a 
lower bound. Knowing the equilibrium value (in terms of the various parameters) for R* n 
in the non-disease model gives another way of calculating d n as the value is constrained 
once a n , a c , /i n , /i c are known. 
Inverse problem 

An experiment on rhesus monkeys in which intact inner segments were transplanted into 
the retina give us estimates on the renewal of rods and cones |26j . The specific data, 
which also gave us the renewal estimates, states that the height of the OS discs on the 
rods (given over mm 2 area and as a fraction of the normal height) was .3, .33, and .44 at 
the respective times of 7, 14, and 30 days after implant. In the cones, the same fractions 
of the overall height of the OS were .31, .35, and .48. At the measurement of 150 days, 
the height of the implanted photoreceptors was indistinguishable from that of the normal 
photoreceptors. Assuming that the observations on this scale hold across the retina and 
having known estimates of the numbers of rods and cones in the retina, we have 8 data 
points that are the averaged over the 12 rhesus monkeys from the experiments. We also 
know the number of RPE at steady state and this gives us a final value for the trophic 
pool. Additionally, the same observations of (26] give observed values of the rod renewal. 

The intervals given represent the range of the four best fit values (using Maximum 
Likelihood Estimator [MLE]) utilizing the average total number of rods and cones in the 
human eye as 91 x 10 6 and 4.5 x 10 6 , respectively [15]. As there is approximately 1 RPE 
cell for every 45 photoreceptors, we estimate the total RPE in the retina at 2.1 x 10 6 . 
Substitution of these quantities gives us the parameter values of our system; see Table [TJ 

The best fit graph is given in Figure HI 
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For the diseased model (pQ), we approximate (3 m = (3 n , d m = d n , a m = a n as first 
approximations since we have not found data on any of those values. To estimate the 
mutation rate, m, we estimated the total number of mutated rods based on the 14 year 
progression of the disease (patient first noticed symptoms at age 10, and died at age 
24, with the postmortem eye analyzed [82] . In doing so, we estimate a mutation rate of 
m = 3.68 x 10- 7 /day. 

Appendix 3: Other Stable Solutions 

As mentioned previously, we have observed stable periodic solutions for certain (unre- 
alistic) parameter values. Analytically, we have shown that no non-degenerate Hopf 
bifurcations (A = iu,u > 0) occur for any of the solutions Ei, ■ ■■ ,E? for any choices 
of parameter values (realistic or unrealistic). For V — 1.5 and the unrealistic values 
a n = 0.0028, a c = 0.00011, a m = 0.0027, m = 0.01, pt n = 0.4, /x c = 0.15, /i m = 0.4, 
d n = 0.0002, d m = 0.00215, k = 0.0001, n = 0.015, (3 m = 0.015, 7 = 0.037, a stable 
limit cycle exists in the full 4-dimensional space; see Figure [5j For the same parameter 
values except for /i c = 0.025, d m = 0.0055, k = 0.001, m = 0.1, we observe a stable limit 
cycle in the 3-dimensional subspace R n = 0; see Figure EJ Following these solutions with 
MatCont has not yet allowed us to determine the bifurcation that gave rise to either of 
these limit cycles. 
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Table 1: Summary of parameters for non-diseased system taken from best 6 MLE results. 
The equilibrium values for rods, cones, and trophic (RPE) are taken as 91 x 10 6 , 4.5 x 10 6 , 
2.3 x 10 6 , respectively. 



Variable 


Description 


Source 


T G [1.42, 1.56] 

k = 5 x i(r 7 

Pn = 1.0 X 10" 9 

7 G [4.55, 5.44] x 1(T 8 
d n G [0.90, 1.40] x 10- 11 

^ n 6 [<p> To] 

^ c e b' 9T5] 
a n G [4.3, 4.6] x 10~ 8 

a c G [4.5, 4.8] x 10~ 8 


Total inflow rate into the trophic factors 
Limiting capacity of trophic factors 
Constant per cell rate at which R n withdraws T 
Constant per cell rate at which C withdraws T 
Constant per cell rate at which R n helps C 
OS shedding rate of i? n 
OS shedding rate of C 
Per cell renewal rod rate 
Per cell cone renewal rate 


MLE, [26] 
MLE, [26] 
MLE, [26] 
[26] 
[IS [26] 

[m [22] 
[nana he] 

MLE, [2B] 
MLE, 
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Figure 1: Flow- diagram of photoreceptor interactions affecting RP. Here we assume two 
distinct classes of rod cells depending on whether or not they have had their functionality 
compromised in any way as a result of the genetic mutation and distinguish them as 
normal and mutated rods, R n and R m respectively. For these different types of rods, 
the rate at which they draw from and utilize the trophic factors and their contribution 
of Rod-derived Cone Viability Factor (RdCVF) may be different HSJ ESI EE2 EH] • In this 
work, we do not attempt to distinguish between the causal mutations of RP. T represents 
the RPE cells, which mediate the passage of growth factors, metabolites, ions, etc. to the 
photoreceptors. Fractions of R n , R m , and C (cones) are interpreted as the outer segment 
(OS) discs whereas fractions of T are viewed as the nutrients, etc. that are brought to the 
photoreceptors via the RPE|28l 138] . 
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D =2.45x1 6 D =2.34x1 6 




2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 

D T xio 6 D T X10 6 

Figure 2: Possible paths to blindness in two hypothetical patients with cone- 
rod RP. The planes illustrate the regions where various stages of RP are present in two 
identical patients (i.e., the same D n , 7, etc.) except for certain aspects of their cone 
shedding to renewal ratio, D c . Both D c values (as well as all other parameter values) are 
within the range of realistic values. The blue asterisk represents the parameter values from 
the best fit (lowest MLE). The solid dark lines that separate the regions of parameter space 
correspond to transcritical bifurcations of the Sample paths 1 and 2 in both patients 
trace the same two sets of values for D m and D?- Path 1 in patient 1 (D c = 2.45 x 10 6 
and is along D rn « 2.45 x 10 6 ) illustrates the path E7 (all photoreceptors present) — > E 6 
(cones die off) — >• E 2 (all rods die off resulting in complete blindness) that corresponds to 
a cone-rod RP path. The same sample path in patient 2 D c = 2.34 x 10 6 is represented 
by Ei — > Ei and corresponds to simultaneous rod-cone death. (When varying only two 
parameters there is a very narrow region of Eq stability, suggesting that C dies out first; 
when varying three parameters, we can obtain E 7 — > E 2 .) The steady state solutions E^ 
undergo transcritical bifurcations and change continuously throughout the plane so that 
near a bifurcation set (represented by the bold lines/curves) two steady state solutions 
are "close" in value. For example as we follow the path E-j — > E 2 the cones and both 
types of rods decrease continuously as the parameter changes and just before crossing the 
bifurcation line/curve into region E 2 all photoreceptors are nearly zero, corresponding to 
a patient being almost blind. Sample path 2 corresponds to the same progression in both 
patients; E 7 (all photoreceptors present) — > E 5 (normal rods die off ) — > E3 (cones die 
off) — > E 2 (mutated rods die off giving total blindness) but the range of parameter values 
in each stage (assuming identical reduction in trophic factors and in shedding to renewal 
ratios) is different. 
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Figure 3: Possible paths to blindness in a patient with rod-cone RP. In contrast 
to Figure [2} the left picture fixes the ratio D m and instead varies D c and Dt (with D n 
the same in both pictures and all parameter values within their realistic ranges). The 
solid dark lines that separate the regions of parameter space correspond to transcritical 
bifurcations of the E^. A particular progress of the common rod-cone RP is illustrated 
by the Sample Path. However, a more likely scenario is that three (or more) parameters 
might change in order to give the common sequence E 7 — >• E4 — >■ E 2 . In the right picture, 
we decrease the three ratios D m , D c , and Dt to obtain the same clinically observed 
sequence of stages of RP. (The asterisk is at the intersection point of each plane and the 
numbers in the respective plane are the subscript of the stable steady state solution Ei.) 
At any given point, it will be the case that multiple options exist to slow the progress of 
RP by changing a combination of D m , D c , Dj> and thereby altering a given path in this 
parameter space. 
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Table 2: Summary of parameters for system ([I]). 



Variable 


Description 


Units 


r 


Total inflow rate into the trophic factors 


1/t 


k 


Limiting capacity of trophic factors 


1/Tt 




Constant per cell rate at which R n withdraws T 


l/R n t 


Pm 


Constant per cell rate at which R m withdraws T 


1/R m t 


7 


Constant per cell rate at which C withdraws T 


1/Ct 


d n 


Constant per cell rate at which R n helps C 


l/Ct 


dm 


Constant per cell rate at which i? m helps C 


l/Ct 


Qi 


Per cell conversion factor 


Rn/T 


Q2 


Per cell conversion factor 


Rm/T 


93 


Per cell conversion factor 


C/T 


A*n 


Per cell death rate of R n 


1/t 


/^m 


Per cell death rate of R n 


1/t 


He 


Per cell death rate of C 


1/t 


m 


Per cell mutation rate of R n 


1/t 




Per cell renewal rate /3 n qi 


1/Tt 




Per cell renewal rate (3 m q2 


1/Tt 


a c 


Per cell renewal rate 751 


1/Tt 
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x 10 7 Best Fit Model Using Maximum Liklihood Estimator (MLE) 
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Figure 4: Best least squares fit (MLE) for inverse problem using data from 

[26] . The quantities in the trophic equation are unknown and all but one of the other 
parameters have a range of values from the experimental literature. The five circles in 
each graph are the data points at 7, 14, 30, 150, and 250 days, with only the last two 
points known in the RPE. The specific values used are a n = 4.5 x 10~ 8 , a c = 4.7 x 10 -8 , T = 
1.48, k = 5 x 10" 7 , 7 = 4.88 x 10" 8 , pi n = 1/9.5, /i c = 1/9, /3 n = 1 x 10~ 9 , d n = 1.29 x 10" 11 . 
The assumed steady state values of R* n w 91 x 10 6 , C* « 4.5 x 10 6 , T* w 2.3 x 10 6 for the 
normal rod OS, cones OS, and RPE, respectively. 
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Table 3: Steady states and corresponding stages of RP. Each of the steady state solutions 
represents a (stable) stage of RP that can be attained for some set of parameter values 
with the exception of E\, which is never stable. Many different sub-types of RP can be 
identified depending on the stable equilibrium that occur as parameters are changed. 



Hi 1 Hi i 

Steady State 


Stage of RP 


RP sub-type 


Physiological Description 


E 7 


Early, Mid-late 


All types 


AT 1 J iii j 

Normal rods, mutated rods, 
cones, trophic factors all exist 


E& 


Mid-late 


Cone- Rod 


Cones have died off; both types 
of rods, trophic factors remain 


E 5 


Mid 


Rod- Cone 


No normal rods; mutated rods, 
cones, trophic factors remain 


E 4 


Mid-late 


Rod- Cone 


all rods have died; cones, 
trophic factors remain 


E3 


Mid-late 


Cone- Rod 


Cones and normal rods have died; 
mutated rods and trophic factors remain 


E 2 


Late 


All types 


Complete blindness with no rods or 
cones; trophic factors remain 


E l 


NA 


NA 


No rods, cones, or trophic factors 



Plot of stable limit cycle with R ,R ,C,T>0 Plot of stable limit cycle with R ,R ,C,T>0 
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Figure 5: Stable limit cycle observed in the full 4-dimensional space. The 

parameter values at which this stable limit cycle exists are V = 1.5, a n = 0.0028, a c = 
0.00011, a m = 0.0027, m = 0.01, fi n = 0.4, \i c = 0.15, /i m = 0.4, d n = 0.0002, d m = 
0.00215, k = 0.0001, /3 n = 0.015, (3 m = 0.015, 7 = 0.037. Note that the axes are not in 
millions but are instead the number of OS and RPE (original definition). 
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Plot of stable limit cycle with R =0, R ,C,T>0 



Plot of stable limit cycle with R =0, R ,C,T>0 




Figure 6: Stable limit cycle observed in the 3-dimensional subspace, R n = 0. 

The parameter values at which this stable limit cycle exists are T = 1.5, a n = 0.0028, 
a c = 0.00011, a m = 0.0027, m = 0.1, /i n = 0.4, fi e = 0.025, /i m = 0.4, d n = 0.0002, 
d m = 0.0055, k = 0.001, (3 n = 0.015, (3 m = 0.015, 7 = 0.037. Note that the axes are not 
in millions but are instead the number of OS and RPE (original definition). 
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